Modified Fowler-Nordheim Field-Emission Formulae from a Nonplanar-Emitter Model 
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Field emission formulae, current-voltage characteristics and energy distribution of emitted elec- 
trons, are derived analytically for a nonplanar (hyperboloidal) metallic emitter model. The tradi- 
tional Fowler-Nordheim formulae, which are derived from a planar emitter model, are modified, and 
the assumption of the planar emitter in the F-N model is reconsidered. Our analytical calculation 
also reveals the backgrounds of the previous numerical discussion by He et al. on the effect of the 
geometry of emitter on field emission. The new formulae contain a parameter which characterizes 
the sharpness of the hyperboloidal emitter, and experimental data of field emissions from clean 
tungsten emitters and nanotip emitters are analyzed by making use of this feature. 

PACS numbers: 00.00.- 
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I. INTRODUCTION 

The Fowler-Nordheim (F-N) theory 1 ! 2 ; 3 ' 4 is one of the 
most important theories of electron field emission. It 
describes experiments, i.e., current -voltage characteris- 
tics of field emission currenljiSii and energy distribu- 
tion of field emitted electrons^ quite well. It should be 
noted, however, that the surface of emitter is assumed 
to be planar in their modeU*2ii 4 although actual emit- 
ters are not planar literally, and the famous F-N formula 
I/V 2 oc exp(— A/V), for example, is derived under this 
assumption. 

This assumption may be considered to be justified 
since the emission area on the surface of emitter may be 
limited to a so small region at the apex of the emitter that 
the area can be regarded as planar. There is however no 
rigorous verification of this expectation. Furthermore, 
emitters are becoming sharper and sharper nowadays, 
e.g., nanotip^^ nanotube^ etc., for which the planar 
emitter model is inappropriate. It is hence worthwhile 
to see how the field emission formulae are modified when 
the (three dimensional) geometry of emitter is taken into 
account. 

Attempts to incorporate geometrical effects into 
the F-N theory have already been made by several 
groupsi 8 i 9 i 10 i 1:L He et al. derived bias fields and image 
charge potentials exactly for the tips shaped like cone, 
hyperboloid, paraboloid, and sphere on cone^ and dis- 
cussed the effect of the geometry of the emitters on field 
emission: 9 They numerically obtained (i) field emission 
currents from the apex of those emitters and (ii) (nor- 
mal) energy distributions of emitted electrons, and they 
concluded that (i) the current-voltage characteristics for 
the nonplanar (and rather sharp) emitters can be fit- 
ted with the relationship I/V 2 oc exp(-A/V - B/V 2 ), 
which is different from the one by the F-N model, I/V 2 cx 
cxp(— A/V), and that (ii) the energy distributions for the 
nonplanar models are much wider than that for the F-N 
planar model. 



In this article, we successfully derive analytical formu- 
lae for current-voltage characteristics and energy distri- 
bution for a hyperboloidal emitter model (although the 
image charge effect is neglected). Even though our for- 
mulae are valid only for conventional emitters whose tip 
radii of curvature are of the order of 100 nm, they still 
reflect the geometry of the emitter and the traditional 
F-N formulae are modified. Furthermore, our analyti- 
cal calculation enables us to discuss the assumption of 
the planar emitter in the F-N model, and also helps us 
understand the backgrounds of the above-mentioned nu- 
merical discussion by He et al. And at the final part 
of this article, experimental data of field emissions from 
nanotip emitters^ are analyzed by making use of a pa- 
rameter, which characterizes the sharpness of the hyper- 
boloidal emitter and never comes in the ordinary F-N 
theory based on the planar emitter model. 



II. HYPERBOLOIDAL EMITTER MODEL 

Our model is illustrated in Fig. ^ It ls composed of 
a hyperboloidal metallic emitter and a planar collector. 
The work function of the emitter, d>, is assumed to be 



Emitter (Work Function: (j)) 

\ Distance: L 



- Collector 



Voltage: Vq 



FIG. 1: A model composed of a hyperboloidal metallic emitter 
and a planar collector. A typical set of parameters in an 
actual experiment is L ~ 5 cm, Vo ~ 3 kV, and <j> ~ 4.4 eV 
(tungsten). The work function of the emitter, <j>, is assumed 
to be uniform. 
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uniform on the surface of the emitter. 

It is implicitly assumed in the ordinary F-N the- 
oj-jjliii that the emission area on the surface of the 
emitter, from which electrons are emitted, is limited to 
a very small region at the apex of the emitter and can 
be regarded as planar. The surface of the emitter is thus 
modelled by a plane, and tunnelling of electrons through 
a one-dimensional potential barrier at the surface is dis- 
cussed. In this article, on the contrary, we duly take the 
shape of the emitter into account and proceed without 
such an assumption. The size of the emission area is 
also what is to be clarified as a result of the following 
calculation. 

We work in an orthogonal curvilinear coordinate sys- 
tem, i.e., an ellipsoidal coordinate system (u,v,p) [Fig. 
121 a)]. It is related to the cylindrical coordinate system 
(p, <p, z) by 



{z = — a cos u cosh 
p = asinusinh?;, 

z + ip — —a cos(u + iv), 



(2.1) 



(2.2) 



where the ranges of the coordinates are < u < tt, 
< v < oo, and < tp < 2tt [Fig. H^b)]. A surface 
v = const, is an ellipsoid, and a surface u = const, is a 
hyperboloid. The expression 1)2. 2JI with an analytic func- 
tion shows that this mapping is a conformal one and the 
ellipsoidal coordinate system in Fig.[2fa) is an orthogonal 
one. 

We attach the ellipsoidal coordinate system in Fig.^a) 
to the model in Fig. such that the plane u = tt/2 is on 
the collector and a hyperboloid u = uq is on the surface 
of the emitter. The parameters a and uq are related to 
the distance from the apex of the emitter to the collector, 
L,by 



L = a cos Mo 

and to the radius of curvature of the tip, R, by 

sin 2 uo T 9 
R = a = Ltan «o. 

COS Uq 



(2.3) 



(2.4) 
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(b) 



FIG. 2: An ellipsoidal coordinate system (u,v,ip). This is a 
3-dim. orthogonal curvilinear coordinate system. The whole 
plane in (a) is mapped onto the region < u < tt in (b) 
through Eq. 12. in for each p. 



TABLE I: Radius of curvature R vs uo for L = 5 cm. 



MO 


0.0010 


0.0015 


0.0020 


R (nm) 


50 


110 


200 



uq is the parameter which characterizes the sharpness of 
the emitter. Typical values of it are shown in Table U It 
should be noted here that this parameter is not contained 
in the ordinary F-N theory since the F-N theory is based 
on the planar emitter model. 

In this coordinate system, the Laplacian V 2 is given 
by 



1 d 







1 d 



d 



V = — — a sin u sinh v — 1 — — — asinitsinhu— — 



y/gdu 



du 
d 2 



Vgdv 



a 2 sin 2 u sinh 2 w dip 2 



dv 
(2.5) 



with yjg = a (sin ucosh v + cos usinh v) sin u sinhw, 
and the Laplace equation \7 2 V = with the bound- 
ary conditions V = on the surface of the emitter and 
V = Vq on the collector plane, where Vb (> 0) is the 
applied voltage, is reduced to 



d ■ d mi \ 
— sin it— Viu) 

du du 



0, 



V(u ) = 0, 
V{tt/2) = Vb. 



(2.6) 



It is then easy to obtain the electric potential V(u) be- 
tween the emitter and the collector: It reads 



V(u) = V 1 - 



In cot(u/2) 
In cot(ito/2) 



(2.7) 



An equipotential surface is a hyperboloid u — const., 
and an electric line of force is a curve along an ellipse 
v, ip = const. The potential energy of an electron in the 
electric field, 



U(u) = -eV(u), 



(2. 



U(z,p = 0) 



U{z,p = 0) 




1 nm 




FIG. 3: Potential energy of an electron, U in Eq. 12.81 
with 12.71 . on the z-axis for L ~ 5 cm, R ~ 100 nm, and 
Vb ~ 3kV. The ?7-axis is drawn at the apex of the emitter. 
The potential near the apex in (a) is enlarged in (b). (a) Due 
to the sharpness of the emitter, a strong field is realized near 
the apex of the emitter, and (b) the thickness of the tunnelling 
barrier is of the order of 1 nm. 
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where — e (< 0) is the charge of the electron, is drawn 
in Fig. |21 for p = 0. In this article, we neglect the image 
charge potential. 



III. GENERALIZATION OF THE F-N THEORY 
FOR THE HYPERBOLOIDAL EMITTER MODEL 

Generalizing the F-N theory>ii2iiii we calculate the 
field emission current I for the hyperboloidal emitter 
model by the equation 

2ir oo oo E 

I = e J dip J dv J dE J dWN(W,E,v)D(W,v). 

— oc — oo 

(3.1) 

E and W are the "total" and "normal" energies of 
an electron, respectively^ N(W, E,v) is called supply 
function^ and N(W, E, v) dWdE dv dtp gives the number 
of electrons incident on the area dv dip around the posi- 
tion (u, tp) on the face of the potential barrier u = uq per 
unit time, with the total energy within the range E to 
E + dE and the normal energy within the range W to 
W + dW. D(W, v) is the probability of an electron pen- 
etrating the potential barrier at (v, p) with the normal 
energy W . We obtain the emission current / by sum- 
ming up the number of emitted electrons per unit time, 
N(W, E, v)D{W, v) dWdE dv dp, for all possible energies 
(W, E) and for all positions (v, p>) on the surface of the 
emitter, i.e., by Eq. (Q) . 

It should be noted in this formulation that the radius 
of curvature of the emitter is assumed to be much larger 
than the de Broglie wave length of an electron (~ 0.1 nm) 
since electrons are treated as localized objects. (For con- 
ventional emitters, the radii of curvature are of the or- 
der of 100 nm.) Under this assumption, the surface of 
the potential barrier u = uq seems planar to an elec- 
tron, and the other surface of the barrier (the end of the 
tunnelling region) u — ut, which is separated from the 
surface u = uq only by the distance of the order of 1 nm, 
also seems planar and parallel to the surface u — uq. This 
observation leads us to the following explicit formulae for 
D(W,v) and N(W,E,v). 

For the barrier penetration probability D(W, v), let us 
consider the Schrodinger equation 



approximation^ to 



V 2 + U(u) 1>(u, v, <p) = Ei>(u, v, ip) (3.2) 

2m 



a sin u sinh v 



dOY 
du ) 



V9 

-- 2m[E-U(u)} 



dO 

dv 



a 2 sin 2 u sinh 2 v 
(3.3) 



It is expected in the tunnelling region that there does 
not exist oscillating-wave mode in the u-direction nor 
damping-wave mode in the v-direction (along an equipo- 
tential curve). It might be hence reasonable to assume 
that 



d9 R 
du 



0, 



80! 
dv 



0, 



(3.4a) 



and 



a sin u sinh v ( dOi 



V9 



a sin u sinhw / dO 



du 



2m[U(u)-W], (3.4b) 



2m(E-W), 



yfg \ dv J a 2 sin 2 Msinh 2 w 

(3.4c) 

where W (< E < U) is regarded as the normal energy. 
We obtain from Eqs. I|3.4b|) and 13.4c|) 



/ ds y/2m[U(u)-W], 



(3.5a) 



9 R (u,v)~ d£\ 2m(E-W) 



1>Z 



a 2 sin 2 u sinh 2 w 

(3.5b) 

ds and dl are the length elements along the curve v, 
p = const, and the curve u, tp = const., respectively, 
given byi£ 

ds = <z\/ sin 2 u cosh 2 u + cos 2 it sinh 2 u du, (3.6a) 



d£ = aV sin 2 u cosh 2 u + cos 2 m sinh 2 w dv. (3.6b) 

One can confirm the consistency of the solution l|3.5|) 
with Eq. (|3.4al) under the situation of small curvature 
and short tunnelling length mentioned in the previous 
paragraph. The barrier penetration probability D(W 1 v) 
for an electron at (v, p) is thus given by 



D(W,v)^e- 2S ^ h , 



(3.7a) 



with the Laplacian V 2 given in Eq. (j2.5|l and the po- 
tential barrier U(u) in Eq. (j2.8jl with l|2.7jl . (m is 

the mass of an electron.) Decomposing the wave func- 
tion as ijj(u,v,ip) = e i®{u,v)/h e ip^/h^ where = 

Or(u, v) + i©j(u, v) is a complex- valued function and p v 
a real constant (angular momentum around the z-axis), 
the Schrodinger equation l|3.2|) is reduced in the WKB 



S(W, v)= ds y/2m[U(u)-W]. 



(3.7b) 



Note that ut is determined by the equation U (ut) — W 
0. 
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The supply function at each point on the surface of 
the potential barrier is given by the same formula as 
that for the planar emitter model, which is derived in 
Ref. 0, since the surface of the barrier seems planar to 
electrons around that point. Hence the supply function 
N(W, E, v) dv dip reads 



N(W,E,v) dvdip 



1 



2tt 2 ^ 3 e( £ +«/ fc * T + 1 



d 2 a, (3. 



where d 2 a is the areal element on the surface of the emit- 
ter defined by 



d 2 a = pdipd£ with u = uq. 



(3.9) 



p is given in Eq. (|2.1[) and dl in Eq. (|3.6b|) . 

With these components together, we first calculate the 
(total) energy distribution of emitted electrons 4 



P(E)= dip dv / dWN(W,E,v)D(W,v), (3.10) 



able us to achieve analytical formulae. The first one is the 
linear approximation of the potential V(u) in Eq. I|2.7|l . 
which one may realize from Fig. E^b) : 



V(u) 



V 



sin uq In cot(«o/2) 



(u - U ). 



(3.12) 



This is valid if [ut — uo)cotwo — (— W/eVo) costto x 
lncot(wo/2) <C 1, which is satisfied in conventional ex- 
periments: (ut — uq) cot Mo ~ 10~ 2 for R ~ 100 nm, 

L ~ 5 cm, Vo ~ 3kV, and W 4> 4.4 eV. (In 

this section, we will often demonstrate the validity of ap- 
proximations with this set of parameters. We hereafter 
call it "case I" for short. Of course, this is not the only 
situation to which our formulae are applicable.) In this 
regime, one can evaluate the action S(W,v), defined in 
Eq. (|3.7bl) . for vj sin«o <C 1 (near the apex of the emit- 
ter) as 



and then obtain the emission current / by 

oo 

I = e J dEP(E). 



(3.11) 



Before going on any further, let us summarize the dif- 
ferences between our formulation and that of He et al. in 
Ref. @. (i) We neglect the image charge potential while 
He et al. calculated it exactly^ (ii) We count the elec- 
trons emitted from all over the surface of the emitter 
through the u-integration in Eq. (|3.1|) while He et al. as- 
sumed, as in the ordinary F-N theory, that the emission 
area on the surface of the emitter is small enough and is 
regarded as a plane. These two points will be addressed 
in the following. 



A. Energy distribution P(E) 



S{W, v) 



2meVo 



sinuo In cot(uo/2) 



ds \J ut — u 



S (W) 



1 



2 V sinuo 



1 , 

77 1 — 77 sm u 



sm tio 



(3.13) 



where the factor Sq(W) = S(W, 0) is the action along the 
z-axis given by 



So (WO 



2aJ2m (-W) 3 , , , , 

— v - ; r -sin 2 M lncot(Mo/2 ). 3.14) 

3eVh 



Unfortunately, it is not possible to carry out the inte- 
grations exactly, but a few reasonable approximations en- 



In this evaluation, we have expanded the length element 
ds given in Eq. (|3.(ia|) around u — uq and v — as 




(3.15) 



r 



and neglected 0((ut — Mo)cotwo) terms. The term la (|3.10() if (v/ smuo) 4 /8 -C (f/sinwo) 2 /2 is satisfied 
(v/ sinu ) 4 in Eq. I|3.13|l is negligible in the formu- within the region v < sin u / v / 2S Q (W)/h [which is the 
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region where the barrier penetration probability D(W, v) 
in Eq. (|3~7|) is not zero], i.e., 8S (W)/h > 1. Note here 
that the WKB approximation (|3.7(1 is a semi-classical ap- 
proximation and is valid for large S{W,v)/h. Typical 
value of S (W)/h is S (-<j))/h ~ 7.6 in the case I. The 
areal element cPa in Eq. (|3.9|1 is also approximated by 



d 2 a ~ a 2 v sin uq dv dip 



(3.16) 



for u/ sin itg <C 1, and the integration one should carry 
out is now 



ma 2 sin 2 u 



ttTT 



P(E) 

E oc 

- I dW J dvvexp 





1 



(E+<t>)/k B T + 1 

2s (w){i + l. . 

n I 2 \smu 



(3.17) 



where the action Sq(W) is given in Eq. I|3.14|l . We 
thus reach the energy distribution formula for the hy- 
perboloidal emitter model: 



P{E) 



2m<j)LR 2 



1 



3irfi 3 (L + R) e ( B +«/*" T + l 

' 7 r(-l/3,2S (E)/h), 



(3.18) 



where 



(2S /n)V3- 



So = S (-<l>) = 



2Ly/2m(t) 3 



ZeVoPo 

with a modification factor for field strength, 

COS Uq 



A) = 



sin 2 uo In cot(uo/2) 



(3.19) 



(3.20) 



which comes from the sharpness of the emitter, and 

r(z,p) = J dtt^e-* (3.21) 
p 

is the incomplete gamma function. Since r(z,p) has an 
asymptotic expansion 

r(z, P ) =p z ~ 1 e - p 

r n-i 

1 + E -( z 2 ) • • • ( z - n ) + °(\p\~ N ) 



n=l 



(3.22) 



for large \p\r^ the energy distribution (|3.18|) is further 
approximated by 



P{E) 



2m(f)LR 2 



1 



3nh 3 (L + R) e W» T + l 
1 



(2S /ny 



■ exp 



(2 S »/»)|l-|^ 



(3.23) 



when the action 2Sq/H is large. Note that we had al- 
ready assumed SS^/Ti ^> 1 in Eq. (|3.17|) . Noticing the 
high energy cut-off by the Fermi-Dirac distribution and 
the low energy one by the exponential factor in the expan- 
sion 13.22|) . the action Sq(E) has been expanded around 
the Fermi level E ~ E F = -cj) in Eq. (|3.23|) . The energy 
distributions (|3.18|l and (|3.23|) are derived under the con- 
ditions 



(ut — Mo) COt Mo 

and 



L 



35 



R eV /3 2Ry/2rK( l 



< 1 (3.24a) 



8S Q /h > 1, (3.24b) 
namely, they are valid for the parameters satisfying 



1 < 8S Q /h < 



l&R^/ ? 2Tn~( l 

3h 



(3.25) 



The energy distribution (|3.23|) is plotted in Fig. 0] It 
is clear from the expression (|3.18() and (|3.23(l that the 
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FIG. 4: (a) The energy distributions 13.231 for zero and 
nonzero temperatures. The radius of curvature of the emitter 
is R — 110 nm (no = 0.0015). The peak of the distribution for 
zero temperature is normalized to unity, (b) The energy dis- 
tributions 13.2311 for different radii of curvature, uo = 0.0018 
for R — 160 nm and uo = 0.0012 for R — 70 nm. Temperature 
is T = 300 K. The peak of each distribution is normalized to 
unity. Other parameters for both figures (a) and (b) are those 
of the case I. 



action So as well as the temperature T characterizes the 
energy distribution P(E). Roughly speaking, the tem- 
perature T determines the high energy cut-off of the dis- 
tribution [Fig. I^a)], and the action So the low energy one 
[Fig.QJb)]. For a sharper emitter, the action So is smaller 
[see Eqs. (|3.19|) and l|3.20[l ]. and the energy distribution 
P(E) is wider [Fig. ^b)]. This agrees with and supports 
the numerical calculation by He et al. in Ref. |9J Further- 
more, we have an explicit expression for a measure of the 
width of the distribution, AE: 



9 is evaluated as 



AE = 



irk B T 



sm(3Ttk B TSo/1i4>) 



(3.26) 



which is the deviation of the distribution P(E) defined by 
AE = ^(E 2 ) - (E) 2 where (E n ) = J^dE E n P(E)/ 
dE P(E). Remember the integral dx e x l d j 
( e x/fc B T + i) = TTk B T/ sm( Trk B T /d) for k B T < dM. One 
can see from the formula (|3.26(1 that AE is more sensi- 
tive to the action So than to the temperature T. In fact, 
the factor (3Trk B TSo/h<t>)/ sm(3irk B TSo/fi4>) is almost 1 
(about 1.03 at room temperature T ~ 300 K in the case 
I), and AE is approximated by AE ~ hcj)/3So- 



B. Emission area 

At this point, we can discuss the assumption on the 
emission area in the ordinary F-N theory and the theory 
of He et al.fi where the emission areas are assumed to be 
small enough and are regarded as planes. In our theory, 
the emission area is not given by hand but is determined 
by the barrier penetration probability D(W, v): It is the 
region where D(W,v), which is the exponential factor in 
the ^-integration in Eq. (j3.17() . is not zero and is esti- 
mated as the region v < sin uq/ y/2So/h. This seems at 
first sight to mean that the emission area is smaller for 
a sharper emitter and to support the above assumption. 
This is however not true. One should take the curvature 
of the area into account. The relevant measure is not the 
size itself but the solid angle of the area seen from the 
center of curvature, i.e., the angle 29 in Fig. [S] The angle 




FIG. 5: Emission area measured by the angle 28 seen from 
the center of curvature. 



tan 6* 



l Mo COS Mo sinh(sinM / v /25o/?i) ^h/2S 



1 - cos 2 m cosh(sin u / v / 2S ~Jh) 1 - h/AS 

(3.27) 

(for small uq), which shows that, for a sharper emitter, 
the action So is smaller, and the angle 28 is larger. In 
this sense, the emission area for a sharp emitter is not 
planar. Even in a conventional situation as the case I, 
the action is So/Ti ~ 7.6 and the angle is 29 ~ 30°, which 
is not so small that the emission area can be regarded as 
planar. 



C. Field emission current / 

It is also possible to execute the final integration (|3.11|) 
with the expression (|3.18|) if temperature T = 0. The 
emission current / at zero temperature is then 



2me<j) 2 LR 2 



37Th 3 (L + R) (2SbA) 4/3 

x (V(l/3,2SoA) - (2Sb/?i) 2/3 r(-l/3,2V^ 

(3.28) 

If the action 2So/^ is large, however, a current formula 
for finite temperature is available: Integrating the energy 
distribution (|3.23ll . one obtains 



I = 



4me(j> 2 LR 2 ^k B TSo/H 



,-2S„/h 



9irh 3 (L + R) sm{3irk B TS /h4>) {2So/hf 

(3.29) 

These current formulae are again valid under the condi- 
tion H3.25f) . It has been known experimentally and the 
F-N theory has explained successfully that field emission 
current does not depend on temperature significantly^^ 
This is also the case with our hyperboloidal emitter 
model. As already mentioned in the previous subsection, 
the factor 3irk B TSo/h(f> found in Eq. (|3.29(l is small and 
the temperature correction to the current is only about 
3% at room temperature T ~ 300 K. The temperature 
dependence is hence negligible in the emission current I. 

Since the action Sq is proportional to Vq -1 [see Eq. 
(|3.19|) ]. the emission current (|3.29|) offers the follow- 
ing current-voltage characteristics for the hyperboloidal 
emitter: 



-3 ocexp 



A = 



AL^2m(t) 3 
3he/3o 



A 
Vo 



a COSUo ,, anw 

sm i*o In cot (uo/ 2) 



(3.30a) 



Remember the corresponding relationship in the F-N 
theory, I /V a 2 cx exp(— A/Vo)^^ Taking the geome- 
try of emitter into account, we have obtained different 
exponent of Vq on the left-hand side. This is one of the 
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main results of this article. The new exponent is a direct 
consequence of the Gaussian integral over the variable 
V in Eq. H3.17[l . or roughly speaking, of the finiteness 
of the emission area, which comes from the finiteness of 
the radius of curvature of the emitter. It is easily ex- 
pected that the similar situations arise for various shapes 
of emitter, where the integration over a finite emission 
area yields some exponent [not necessarily the same as 
in the formula l|3.30[l but] different from the conventional 
one. The geometrical effect, in general, manifests itself 
in the exponent^ 

It is interesting, on the other hand, that the expo- 
nential factor on the right-hand side of Eq. (|3.30a() is 
the same as that of the conventional one. This, in some 
sense, supports the validity of the planar model, since the 
current-voltage characteristics are mainly dominated by 
this exponential factor but the power on the left-hand 
side is relatively less important. It is not easy to observe 
the difference in the exponent of Vq experimentally^ and 
the conventional F-N theory works well fortunately. 

Our formula (|3.3()|1 however does not coincide with 
the one that He et al. concluded in Ref. @ for the hy- 
perboloidal emitter model, i.e., I/Vq oc exp(— A/Vq — 
B/Vq). What is the reason for this discrepancy or what 
is the origin of the term —B/Vq? The discrepancy seems 
at first sight due to the fact that the image charge effect is 
neglected in our formulation while He et al. calculated it 
exactly. This is not the case, however. In general, the im- 
age charge effect becomes more prominent as the applied 
voltage increases, and the emission current is accordingly 
enhanced. But such an effect is not observed in their cal- 
culation. The additional —B/Vq term in the theory of 
He et al. comes not from the image charge effect but from 
the bias field potential. In fact, if one retains the term 
of the order of (ut — u )cotuo — (L/R)(— W/eVoflo) m 
the action integral (|3.14l) . which we have neglected there, 
such a correction appears: 



S (W) 



2Ly/2m(-W) 3 



a- 



L (-W) 



(3.31) 

where a is a positive constant of the order of 1^ This 
correction comes in when the emitter is so sharp that 
the linear approximation of the bias field potential V(u) 
in Eq. (|3.12(l is not valid. Note the i?-dependence of 
(ut - u )cotu ~ 3S Q (W)/2Ry / 2m{-W). It should be 
noted, however, that emitters with very sharp radii do 
not fit in well with our formulation as mentioned in the 
second paragraph of this section. 



IV. ESTIMATION OF EXPERIMENTS 



experimental data to be fitted in the following were ob- 
tained for tungsten emitters in the experiments whose 
details are presented in Ref. 0- In those experiments, the 
parameter L is about L ~ 5 cm. 



A. Current— voltage characteristics 

Experimental data of current-voltage characteristics 
for different emitters A-D are shown in Fig. El Not only 
those for normal clean emitters but also those for nan- 
otip emitters^ which are fabricated on the apexes of the 
formers are plotted in the same figures. The vertical axis 
I/Vq in Fig. is different from that of the ordinary F- 
N plot, I/Vq, and each series of data is well fitted by a 
straight line, i.e., by the formula H3.30J1 . 

As already mentioned in the previous section, fittings 
of the data with the conventional F-N formula also work 
well similarly to the ones with our formula shown in 
Fig. El Our formula (|3.3()(l however contains the parame- 
ter ito, which characterizes the sharpness of the emitter, 
while the conventional one does not. We can hence es- 
timate the radii of curvature R of the emitters by the 
fittings, which are listed in Table ITT1 It should be noted 



> 



10" 



10" 



10" 




2 4 6 

1/Vo (lO^V" 1 ) 



> 



10" 



io- 



10" 




1/V (lO^V" 1 ) 



In our formulae, we have the parameter uq, which char- 
acterizes the sharpness of the hyperboloidal emitter and 
is not contained in the ordinary F-N theory based on 
the planar emitter model. Finally, let us try to analyze 
actual experiments by making use of this feature. The 



FIG. 6: Experimental data of current-voltage characteristics 
for four different pairs (A-D) of emitters: Dashed lines are for 
normal clean emitters, and solid lines for nanotips fabricated 
on top of the clean ones. 6 Each series of data is fitted by the 
formula I3.3UH . 
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TABLE II: Radii of curvature of the emitters estimated by 
the fittings in Fig. HJ with (j> = 4.4 eV and L = 5 cm. 





Clean 


Radius of curvature R (nm) 
Nanotip 


A 


154 


99 


B 


86 


43 


C 


201 


125 


D 


108 


58 




(a) (b) 

FIG. 7: (a) Radius of curvature of a clean emitter and (b) ef- 
fective radius of curvature of a nanotip emitter. 

there that the radii of curvature for the nanotips are not 
exactly those of the tips of the emitters: They are just 
the radii of curvature of effective clean emitters which 
emit the same currents as the actual emitters with nan- 
otips (Fig. 0| . One can see from this estimation, however, 
that a nanotip fabricated on top of a clean emitter has 
such an effect that it reduces the radius of curvature of 
the emitter to the extent of a half of the original clean 
one. 



B. Energy distribution 

Energy distributions for normal clean emitters are also 
fitted well by the formula l|3.18[) and (|3.23l) as shown in 
Fig. [SJ The parameter uo, i.e., the radius of curvature 
R of emitter, is again obtained by the fitting. The ra- 
dius of curvature of the clean emitter A in Fig. El which 
is estimated with the data for different applied voltages 
Vq , is presented in Table Ifffl Although the radius of 
curvature R should be independent of the applied volt- 
age Vq, of course, the estimated value decreases slightly 
with Vq. This is due to the image charge effect. This 
effect is neglected in our formulation, and the enhance- 
ment of the emission current due to it is renormalized into 
the parameter Uq, i.e., into the sharpness of the emitter. 



TABLE III: The radius of curvature R of the clean emitter A 
in Fig.|S|estimated by fitting energy distributions for different 
applied voltages Vb with the formula 13.2311 . The parameters 
are fixed <j> = 4.4 eV, T = 300 K, and L — 5 cm. 



Vo (kV) 


3.8 


4.0 


4.2 


4.4 


R (nm) 


114 


113 


111 


95 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 



E-E F (eV) 




-2.0 -1.5 -1.0 -0.5 0.0 
E-E F (eV) 



FIG. 8: Experimental data (dots) of the energy distribution 
for the clean emitter A in Fig.|S|at Vb = 3.8 kV, fitted (lines) 
by the formula l)3.23|l in the log-linear plot (b) with the pa- 
rameters Ep = — <j> = -4.4 eV and T = 300 K fixed. The 
linear one (a) is drawn with the same parameters as in (b). 

The enhancement is a slowly increasing function of Vb 
as already mentioned in the previous section, and the 
estimated radius of curvature of the emitter accordingly 
decreases with Vb- One hence has to be careful about 
the fact that the image charge effect is included in the 
parameter uq. 

V. CONCLUSION 

In this article, we have derived field emission formu- 
lae from a hyperboloidal emitter model, i.e., current- 
voltage characteristics (|3.28[) . (|3.29|) . and (|3.3U[1 . and en- 
ergy distribution of emitted electrons (|3.18|) and i|3.23[l , 
which are valid under the condition 13.25fl . Reflecting 
the geometry of the emitter, the traditional F-N for- 
mulae, which are derived based on the planar emitter 
model, are modified. The current-voltage characteris- 
tics in the F-N theory, I/V 2 cx exp(— A/V), is replaced 
with I/V 3 ex. exp(— A/V), for example. The geometri- 
cal effect manifests itself in the exponent of Vq on the 
left-hand side. 

We have also addressed and reconsidered the assump- 
tion of the planar emitter in the F-N theory. An esti- 
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mation of the spread of emission area based on the for- 
mula l|3.27[) shows that the area cannot be regarded as 
planar even for a conventional emitter. 

Furthermore, our analytical calculation has revealed 
the backgrounds of the conclusions drawn by He et al.^ 
which are based on a numerical calculation: The origin 
of the correction term in their current -voltage character- 
istics has been clarified, and the dependence of the width 
of energy distribution on radius of curvature of emitter 
has been explained. The concise formula H3.26[l for the 
width of energy distribution might be useful in practical 
experiments. 

And finally, we have attempted to analyze experimen- 
tal data of nanotip emitters^ by making use of the param- 
eter Mo, which characterizes the sharpness of the emitters, 
and clarified an effect of a nanotip fabricated on top of 
a normal clean emitter: The effective radius of curvature 
of the emitter is considerably reduced. 



One should note, however, that our formulae do not 
explain the characteristics peculiar to nanotip emit- 
ters exactly. Actually, energy distributions for nanotip 
emitters)** for example, cannot be fitted by our for- 
mula H3.18fl an d (|3.23l) . One of the reasons for this 
is that electrons cannot be treated as localized objects 
in nanotips, whose radii of curvature are of the order 
of 1 nm. For a rigorous description of field emission 
from today's nanoscale tips, a fully quantum treatment 
is required 
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